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Abstract. Triple donor devices have the potential to exhibit adiabatic tunneling 
via the CTAP (Coherent Tunneling Adiabatic Passage) protocol which is a 
candidate transport mechanism for scalable quantum computing. We examine 
theoretically the statistics of dopant placement using counted ion implantation 
by employing an analytical treatment of CTAP transport properties under 
hydrogenic assumptions. We determine theoretical device yields for proof of 
concept devices for different implant energies. In particular, wc determine a 
significant theoretical device yield (~ 80%) for 14 keV phosphorus in silicon with 
nominal 20 nm spacing. 



1. Introduction 

Quantum computer (QC) architectures compatible with conventional silicon 
processing technologies would seem to be at an advantage over other schemes. This 
advantage is due to the perceived ability to leverage conventional silicon processing 
techniques. Of particular interest is the Kane solid-state quantum computer [1] which 
is based around nuclear spins of phosphorus ( 31 P) donors in an isotopically pure 
28 Si matrix (Si:P). As all of the controls and couplings are essentially electron spin 
mediated or effected, an electron spin version version was proposed by Hill et al. [2] 
that took advantage of global control properties to access the faster electron, rather 
than nuclear, spin gate times. 

The Si:P QC still faces substantial technical challenges before it can become 
a reality and effective means of transport such as that proposed by CTAP 
(Coherent Tunneling Adiabatic Passage) [3] offer significant flexibility in overcoming 
these challenges. Amongst these challenges are the expected penalty for linear 
nearest neighbour (LNN) architectures in decreased threshold for fault-tolerant error 
correction due to extra SWAP operations [HE], and with Si:P, valley degeneracy leads 
to sensitivity of most critical parameters to variations in donor positions at the atomic 
level [6j [7] . Both these problems can be ameliorated by the use of effective long-range 
transport which allows the LNN bottleneck to be avoided, and the incorporation of 
defect-tolerant design methodologies. An architecture based around transport via 
CTAP was recently proposed [8] and a threshold analysis performed for a bilinear 
geometry [4]. The use of CTAP transport rails as designed in [8] also avoids a 
complex classical control issue, namely the prohibitively high gate density of the 
original Kane scheme [9|. Although the proof-of-concept structures that we discuss 
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here do not provide the advantages of CTAP rails, demonstrating three-dopant devices 
is a critical step towards eventual scalable structures. The preferred method for 
fabricating longer CTAP chains is hydrogen resist lithography that allows for near- 
atomic dopant placement [10] . This bottom- up process has recently had remarkable 
success in creating nanoscale surface structures [IT] [12] for quantum devices. However 
the demonstration of the CTAP mechanism requires the construction of single dopant 
arrays deep within a silicon substrate. Currently, the quickest route to the construction 
of a proof-of-principle device is the top-down method of single ion implantation |13j as 
it can be adapted to such fabrication with existing technologies. Here we investigate 
this practical method for engineering a three atom device. We show estimates for the 
times required for high-fidelity CTAP using hydrogenic approximations to the tunnel 
matrix elements [14] and ion implanted donor positions calculated using the SFJM [15] 
package. Analytic solutions to the expected time for CTAP based on the adiabaticity 
criterion allow some statistics of expected device yields to be made, which are of 
benefit in guiding experimental investigations. We also discuss some of the limitations 
of our approximations, and point to some of the challenges in identifying complete 
adiabatic pathways to realise the CTAP protocol. 

2. Coherent Tunneling Adiabatic Passage and Adiabaticity 

CTAP is a protocol for the spatial transport of a particle between two-points on a 
quantum chain. In its simplest case this is a three-site protocol with the central chain 
being a single site. CTAP is the direct spatial analog of the well-known STIRAP 
protocol from quantum optics [Hj. CTAP is distinguished from STIRAP in that 
with CTAP, all variations in couplings are effected by direct modulation of the wave- 
function overlaps by either surface gate control or well-proximity. With STIRAP the 
tunnel matrix elements arc- strictly electromagnetic field driven 17 . In addition to 
the electronic transfer exploited for Si:P systems, CTAP has been proposed in atomic 
lattices [18] , Cooper-Pair boxes [19] , quantum-dots [3] [20] , Bose- Einstein condensates 
[21] [22] , spin chains [23] and for photonic transport through coupled waveguides [24] . 
The latter case has also been demonstrated in beautiful experiments by Longhi and 
co-workers [25] [26] which adds significant impetus to the push to demonstrate CTAP 
in systems of massive particles. 

CTAP is also an example of a growing search for direct analogues of existing 
quantum optical effects in quantum electronic systems, and as the degree of coherence 
that can be detected increases one should expect more sophisticated effects to be 
observed. Other than CTAP, there are also proposals for Autler-Townes measurements 
[27l [28] , and coherent-population trapping- like dark state transport protocols [29] [30] . 
One attraction of quantum electronics over conventional quantum optics lies in the 
ability to tailor the Hilbert space by construction, rather than being limited by the 
structures of atomic systems. This flexibility allows new vistas and extensions to be 
explored, for example the multiple-recipient adiabatic passage of Refs. [3TJ[32] and the 
direct spatial analogue of the tripod atom [33] . 

To effect the CTAP pathway requires the adiabatic transformation of one of the 
eigenstates of the system from a known start state to a desired final state. In the 
spatial setting we restrict the control parameters to the tunnel matrix elements which 
are controlled by surface electrodes (gates), analogously to the way in which electro- 
magnetic field intensities are varied in STIRAP. A schematic of the surface structure 
that we are assuming is shown in Figure Qua)] viewed from the top down. The device 
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SIDE VIEW 



Figure 1: (a) Triple dopant, one charge system with surface gate control for realising 
CTAP. The symmetry gates Si and S3 maintain the degeneracy of the end of chain 
states, whilst the variations in the tunnel matrix elements is effected by the barrier 
gates B12 and i?23- (b) top veiw (c) side view SRIM simulations showing the spatial 
probability distribution of 100,000 14keV phosphorus ions implanted into silicon 
through three apertures 10 nm in diameter spaced 20 nm apart. 



structure shown in Figure ( J^a)| is yet to be made and it is the purpose of this paper 
to investigate potential fabrication methods based on presently available technologies. 
Devices of comparable complexity have been fabricated using the counted-ion implan- 
tation techniques described in Ref. [13j . and we also note the demonstration of triple 
dots in GaAs 2DEG structures [3H [351 UH] and gated carbon nanotubes [37J . 

Ideally we would like three donor atoms spaced 20-30 nm apart, 20 nm below the 
surface. The SRIM data in Figures D Pfo) | and [ ^ c) | shows the probability distribution of 
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100,000 14kcV P + ions implanted sucessively through three circular apertures 10 nm 
in diameter. A 14keV phosphorous ion will travel around 20 nm into the substrate 
before stopping but can straggle away from this median position by up to llnm as 
it undergoes collisions with the Si lattice. A less energetic 7keV ion will penetrate 
14 nm below the surface and has 40% less straggle than at 14kcV. The straggle im- 
poses constraints on the accuracy of dopant placement. To quantify the consequence 
of these fabrication perturbations on the ideal donor positions seen in Figure Q [^a)| we 
have used SRIM simulations with a simplified treatment of three-donor CTAP. 

We begin by writing down the Hamiltonian of the three donor one electron 
problem in the three-state approximation (i.e. where we only keep the lowest state of 
the electron localised around the donor). In the basis |1), |2), |3) with onsite energies 
Ei = (assuming fully compensated) and tunnel matrix elements O12 (t) and ^23 (i) 5 
we also include the possibility of next nearest neighbour tunneling (i.e. from |1) to 
|3)) however for the discussion that follows this will be assumed to be zero. The 
Hamiltonian is 

H = fi 12 (*)|2)(l| + 23 (i)|3><2| + fi 13 (i)|3)<l| + h.c, (1) 
and with fli^(t) = 0, the eigenvalues are 

n 23 |i)-n 12 |3) 



v/n? 2 + n 



with energies 



w v = n 12 |i)±^ 2 + ^ 3 |2) + n 2 3|3) 
N/2(n? 2 + n^ 3 ) 

E = 0, (4) 

E ± = ± y^Tn^. (5) 

The CTAP protocol can now be understood quite simply. The idea is to remain in the 
state |X>o}) and to vary the tunnel matrix elements so that at time t = 0, the system 
is in the desired initial state, e.g. \T> (t = 0)) = |1), and at time t = t max the system 
is in the desired final state, e.g. \T> (t = t max )) = |3). Note that although this three 
mode description for the tunneling is obviously a simplification, it still captures all of 
the essential physics of the CTAP protocol, a fact confirmed by recent analyses of the 
CTAP in the triple square well case [38] . 

To effect CTAP, there is clearly a large amount of flexibility to choose the pulsing 
scheme for the tunnel matrix elements. In STIRAP protocols, gaussian or gaussian-likc 
pulses are most commonly employed because of the necessity to turn on the excitation 
(laser pulse) before varying it |16| . However in the solid state, where non-zero tunnel 
matrix elements arise solely due to donor proximity, this necessity is not required, 
and so we advocate the use of the pulses that vary between their extrema at t = 
and t — t max - Such pulses are explicitly stated below and illustrated in Figure In 
Refs. mi [321 error function pulses were employed which have some advantages in 
terms of smoothness of evolution and in nonlinear systems avoid certain complications 
due to eigenstate degeneracy at the ends of the protocol. For simplicity, here we choose 
sinusoids with the sole purpose of making the analytical results clearer. 

To determine whether the system remains in the target state, we invoke the 
adiabaticity criterion, and due to symmetry we define (without loss of generality) the 
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Figure 2: (a) Two possible scenarios for implanted triples, (b) Tunnel matrix elements 
as a function of time. The solid lines correspond to the case that Wyi = W23 for the 
sinusoidal variation as defined in the text, whereas the dotted line corresponds to 
the case that W23 = W12/2. (c) Eigenenergies of the states as a function of time 
with the sinusoidal variation, again solid lines are for W12 = W23 and the dashed to 
W23 = W12/2. (d) the value of the adiabaticity parameter A throughout the process. 
Note that it is maximised at t = t max /2 irrespective of the values of W12 and W23, 
although there are minor differences in the value of A. As expected, the process is 
slightly less adiabatic with the smaller values of W. 



adiabaticity parameter to be between |2?p) and The adiabaticity parameter is 



\E + - E «r () 

and for adiabatic evolution we require i « 1. One should be mindful of the fact 
that the adiabaticity does not translate to a direct measure of fidelity, it is rather a 
measure for when the assumption of adiabatic evolution is justified. Choosing as the 
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form for the (gate controlled) tunnel matrix elements 

n 12 (t) = W 12 sin 2 (-^-), (7) 

n 23 (t) = W 23 co S 2 (-^-), (8) 



2/, 



gives 



irW 12 W 23 sin (f^-) 
A = ywy , g , 

V2t max {W? 2 + W% 3 ) 3/2 

Of importance here is to note that the adiabaticity is a maximum when t = t max /2 
irrespective of the relative values of W± 2 and W 23 . This is significant in the design 
of robust sequences and in affording a quick estimate of the required timescale for 
CTAP operation. To illustrate the energies, tunnel matrix elements and adiabacity, 
in Figure [2] we present characteristic plots for the case that W\ 2 = W 23 and when 
W 23 = W 12 /2. 

In a realistic experiment, we will want to set the time for CTAP given a certain 
desired adiabaticity. So it is more important to rearrange © to determine the value 
for i max that keeps the maximum value of A at or below some threshold, which is 

nW 12 W 23 



V2A(W? 2 + Wl 3 f 2 ' 



(10) 



2 3) 

where we have dropped the sine term because we are evaluating i max with respect to 
the maximum value of A. 

Equation (|10p is particularly useful in gaining insight into the practicalities of 
CTAP. Although we are not yet in a position to place a lower bound on i max (in general 
this must come from a detailed understanding of the higher lying molecular and orbital 
states, limitations on the pulsing electronics, and accessible surface gate voltage), we 
can immediately compare our results with realistic upper bounds. The most obvious 
impediment to very long time scales will be the limits placed by decoherence, which 
will set a maximum length of time over which the protocol can be conducted [3l 140"! I4T] . 

Until now we have neglected the next nearest neighbour tunneling, i.e. ^13. 
Formally there is no CTAP pathway for non-zero SI13 as the Hamiltonian fl} has no 
null (or dark) state, however we can place a good bound on whether neglecting Qi 3 
will be valid by comparing the period for oscillation on the |1) — |3) transition directly 
with the total time for CTAP. We introduce J = r2 13 t max and assert that if J <C 1 
then we may ignore the effect of f2 13 in our analysis. This criterion is helpful for proof 
of concept devices, but may not suffice for full QC applications where more rigorous 
error control is required [5]. 



3. Implications for ion implanted devices 

We have treated the Hamiltonian in an ideal case, without regard to the physical 
nature of our system. If we now turn specifically to phosphorous dopants in silicon, 
with realistic gate controls, then issues of the suppression of tunnel matrix elements 
via barrier control, cross-talk between gates [42] and microscopic details of the dopant 
species become paramount [BJ. All such features must be taken into consideration to 
determine the adiabatic pathway, which is a highly non-trivial task which we do not 
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attempt here. Instead we have focussed on hydrogcnic approximations which allow us 
to gain rapid insight into a variety of implant conditions; we leave the solution of the 
control problem to future work, assuming (as in Rcfs. 3,8) that the gate controls are 
able to completely suppress the tunnel matrix element from the maximum (no-field) 
level to zero. 

To determine the bare (unperturbed) tunnel matrix elements, we use the 
hydrogcnic approximation of Opcnov [14] . an approach used ostensibly in the study 
of singly ionized double-donor structures 43J. These results are applicable to 
donor separations predominantly along the [100] direction in silicon where the valley 
degeneracy has less effect [6||44]. To be explicit, we use 



"^•(SM-t- 1 )' (u » 

E* is the effective Hartrec, and dij is the interdonor separation between donor i and 
j. With these constraints we can immediately write down i max and J as 

d 12 +d 2 3 _ 



ira* B di2d 2 3 exp 



AV2E*A[(F l2 f + (F 23 ) 2 } 

71^12^23^13 exp ^- 



3/2 : 



(12) 



* = } * ^ ,3/2 ~ . (!3) 

V2A[(F 12 j 2 + (F 23 ) 2 } 

F, =^e XP (-|-l)=^. (14) 

A rigorous calculation of t max and J would involve full band structure considerations, 
but these formulae give an extremely efficient mechanism for determining these 
important parameters with relatively minimal computational cost, as is required 
to process the large number of dopant positions that can be obtained using SRIM 
calculations. 

To explore the range of expected yields, we performed SRIM simulations for 
14keV and 7keV phosphorus. Single 14keV phosphorous ions are routinely detected 
entering a silicon substrate |13j and reach a nominal depth of 20 nm below the substrate 
surface. The lateral straggle of an implanted ion limits the precision of dopant 
placement and hence the control of donor seperation required to engineer a suitable 
proof of concept device. Therefore 7keV was also chosen for the greater placement 
accuracy it allows and as a guide to future developement of the detection system. A 
typical oxide thickness of 5 nm was used in the 14keV strategy while in the low energy 
strategy we assumed an oxide thickness of an optimal 1.2 nm. The on-chip detection of 
[T3] guarantees that precisely three ions will be implanted, however it does not ensure 
that there will be one ion per aperture. To ensure this capability a 'step-and-repeat' 
strategy based on existing technologies is currently being integrated with the detection 
system. A similar approach is being followed by Schenkel et al. [45l [46] and also by 
Meijer et al. [47] for the creation of nitrogen-vacancy colour centres in diamond. 

The results of performing SRIM simulations for each implant strategy are shown 
in Figure [3] In addition to the raw straggle data obtained from SRIM, we have 
convolved the distributions with apertures 10 nm in diameter, separated linearly by 
20 nm. For each implant triple, we have then calculated i max and J, although as 
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Figure 3: Cumulative distribution functions showing the fraction of implant triples 
with CTAP times at adiabaticity A = 0.01 for different implant strategies. The solid 
line is for 14 keV P implanted through 5 nm of Si02 on Si; the dashed line is for 7 keV 
P through 1.2 nm of SiC>2 on Si, a practical lower limit to conventional fabrication. 
The horizontal line indicates 10 -9 seconds, a realistic limit on the CTAP transfer 
time due to dephasing. From this we may infer that the 14kcV implant strategy 
has an expected yield around 80%, whilst the 7keV strategy has an expected yield 
around 90%. Note that these yields do not take into account the lower bounds on 
controllability, and some fraction of dopant triples will have tunnel matrix elements 
that are too large to be effectively controlled. 

J is found to be substantially less than one for over 95% of the triples we have not 
shown those results here. Figure [3] then shows the cumulative distribution function 
for the i max at an adiabaticity of A = 0.01, which may immediately be interpreted as 
a measure of yield. Upper limits on i max are given by the decoherence limit. As a rule 
of thumb, we wish to be at least an order of magnitude faster than decoherence for a 
significant proof-of-concept signal [3]. To aid comparisons, Figure [3] shows a horizontal 
line corresponding to 1 ns, which is an order of magnitude less than a realistic expected 
decoherence time (10ns). From this we can read off the population of triples which 
satisfy this criterion for each implant strategy. One should be mindful, however, that 
the expected device yield even without fabrication errors will be less than this, due to 
the fraction of donors that are too close (i.e. have too high tunnel matrix elements) to 
be properly controlled, or to have their tunnel matrix elements adequately suppressed. 
This bound will be determined in large part by the breakdown of the oxide barrier 
and is not addressed here. 

From the results in Figure [3] we can make some initial estimates of the prospects 
for ion implantation strategies, and in particular we find that within the limits of our 
models, we can afford some optimism about the ion implantation as a strategy for 
proving the concept of CTAP. The expected theoretical yield for the standard 14kcV 
strategy is predicted to be around 80%, and for the low energy 7keV strategy it is 
expected to be closer to 90%. 

4. Conclusions 

We have performed analysis of implanted donor triples using SRIM to provide 
guidance to experimental efforts to realise the CTAP (Coherent Tunneling Adiabatic 
Passage) protocol. Our results suggest theoretical device yields around 80% within the 
hydrogenic approximation and not including device fabrication errors. To process the 



Strategies for triple- donor devices fabricated by ion implantation. 



large amounts of data generated by SRIM, we developed new analytical results for the 
required time for high-fidelity CTAP, and applied these metrics to various implantation 
strategies and species. Our methods are unable to address all of the requirements for 
CTAP. which include the need to identify the particular adiabatic pathway for a given 
microscopic location of donors, nor do we go beyond the hydrogenic approximation. 
More accurate modelling is required, (e.g. NEM03D (48l[49j) to theoretically describe 
a given system, but this is far too numerically intensive to explore the very large 
number of configurations explored here. On the basis of our results, we predict 
that ion implantation is a sensible strategy to exploit for proof-of-concept CTAP 
devices, although the overall scaling to the longer CTAP chains required for scalable 
QC [8] is expected to be far less favorable and hydrogen resist lithography methods 
for producing such CTAP chains is the preferred method [10] . 
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